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^H ^ Abstract. Using spectroscopic data presented in Magrini et al. 12UU3II . we have 

analyzed with the photoionization code CLOUDY 94.00 (Ferland et al. I1998II 11 



Planetary Nebulae belonging to the spiral galaxy M 33. Central star temperatures 
and nebular parameters have been determined. In particular the chemical abun- 
dances of He/H, O/H, N/H, Ar/H, and S/H have been measured and compared 
j^ , with values obtained via the Ionization Correction Factors (ICFs) method, when 

available. Chemical abundance relationships have been investigated; in particular, 

• i-H . 

^^ ■ a correlation between N/H and N/0 similar to the Galactic one fHenrv [199011 . 

^^ 

jrt ' and a feeble anti-correlation between O/H and N/O have been found. A gradient 

in O/H across the disc of M 33 is indicatively consistent with the one found from 
H II regions in this galaxy (Vilchez et al ll988t . Further studies in the more exter- 
nal parts of M 33 are however needed to ascertain this point. The present result 
shows that oxygen and helium abundances (with lower accuracy also nitrogen, 
argon and sulphur) can be actually estimated from the brightest PNe of a galaxy, 
even if the electron temperature cannot be measured. We also found that the 
oxygen abundance is quite independent of the absolute magnitude of the PN and 
consequently the brightest PNe are representative of the whole PN population. 
This represents an important tool to measure the metallicity of galaxies at the 
time of the formation of PNe progenitors. 
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1. Introduction 

Chemical abundances in extragalactic Planetary Nebulae (PNe) have been derived so 
far only in a small number of galaxies. Inside the Local Group they are: the irregular 
galaxies LMC, SMC (cf. Leisy et al. l^nOUl Jacoby & De Marco EDHl, and NGC 6822 
(Richer & McCali ri995|l . the dwarf spheroidal and spheroidal galaxies Fornax (Maran et 
aL lT^MjI . Sagittarius (cf. Walsh et al. lTWTIl . NGC 185, and NGC 205 (Richer & McCaU 
I1995|l . the elliptical galaxy M 32 (Stasinska et al. 11998(1 and finally the two spiral galaxies 
M 31 and M 33 that have been studied by Jacoby & CiarduUo IJ1999I hereafter JC99) 
and Magrini et al. I)2()()3[ hereafter M03), respectively. There is a substantial lack of 
information about chemical abundances of PNe in dwarf irregular galaxies of the Local 
Group (e.g. Sextans B. Sextans A, Leo A, IC 10), where no PN has been investigated so 
far. Outside the LG, only PNe in the giant elliptical galaxy NGC 5128 (Centaurus A) at 
the distance of 3.5 Mpc were analyzed by Walsh et al. (19995. 

In most of the previous cases, with the exception of NGC 5128, the abundances have 
been obtained by direct determination of the Te, derived from the [O in] AA 4363/5007 
A or, in the worst cases when the A 4363 A is not measured, using its upper limit, thus 
obtaining lower limits for the [O ill] ionic abundance and then for the total O abundance, 
which was calculated from the ICFs (Ionization Correction Factors) procedure (cf. e.g. 
Kingsburgh & Barlow HlMl hereafter KB94). 

In M 31, in addition to the ICFs method, JC99 have modeled the observed PNe using 
the photoionization code CLOUDY 90 (Ferland et al. ll998|l with simplified assumptions: 
blackbody central stars and spherical nebulae with constant density. They have derived 
abundances and central star parameters for 15 PNe of M 31. In eight cases they could 
compare the results from the nebular models built with CLOUDY with abundances based 
on the ICFs method. They found good agreement between abundances derived with 
these two different methods. From this comparison, CLOUDY resulted to be a useful 
tool to measure chemical abundances, particularly when direct electron temperature 
measurements are not available. 

M03 studied 48 emission line objects in M 33 spectroscopically, and recognized, via 
diagnostic diagrams, 26 of them to be bona fide PNe. They presented the observed fluxes 
for 42 of the 48 objects, deriving chemical abundances with the ICFs method for the 
three brightest PNe. 

In the present paper, we re-analyze with CLOUDY 94.00 (Ferland et al. I199X|I 11 
of the 26 PNe presented by M03 using the available spectroscopic data. The three PNe 
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already studied with ICFs by M03 are also included, in an effort to clarify their observed 
low N/0. 

In Section 2 we summarize the observations and data reduction. In Section 3 the 
derivation of nebular and stellar properties is illustrated. In Section 4 we discuss the 
derived quantities and correlations among them. Summary and Conclusions are in Section 
5. 

2. Observations and data reduction 

The observations were taken at the 4.2 WHT telescope (La Palma, Spain) in October, 
2001 with AF2/WYFF0S, with a setting producing a dispersion of 3.0 A per pixel in 
the spectral range 4300 A to 7380 A. A detailed description of spectroscopic data and of 
the reduction procedure are in M03. 
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Fig. 1. Characteristic spectra of PNe in M 33 obtained by M03 with an exposure time 
of 8.7 hrs with the WHT. 
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3. The modelling procedure 

Using a photoionization model such as CLOUDY, it is possible to overcome the lack 
of direct information on the electron temperature in a PN, due to the faintness of the 
relevant forbidden lines: the [O in] A4363 and the [N ii] A5755 emission lines, which are 
hard to detect in extragalactic PNe. On the other hand, the use of CLOUDY requires as 
input the energy distribution of the ionizing radiation of the central star, its luminosity 
and the nebular geometry. The distance to the object must then be known. In our case 
the distance is available adopting for all PNe in M 33 the distance of the galaxy, equal 
to 840±90 kpc (Magrini et al. I2()()()l hereafter MOO). As with the energy distribution 
of central stars, the usual simplification is to assume them to behave as blackbodies 
(BB), while for a better accuracy stellar atmosphere models are needed (cf. Morisset 
I2004|l . In this respect, an illustrative test was carried out by Howard et al. H1997|l . They 
tested the effects on derived abundances and central star temperatures of a variety of 
model atmospheres and blackbodies in a set of Galactic halo PNe. They found that 
nebular line strengths are relatively insensitive to atmospheric model details and that 
blackbody spectra are by and large suitable to represent the central star continua. On 
the other hand Armsdorfer et al. H2U(J2fl checked in the particular case of the spatially 
resolved round Galactic PN NGC 2438, the effects of modelling a PN central star with a 
blackbody instead than with a model atmosphere. They found that the BB assumption 
leads to underestimating the [O in] A5007 line and to overestimating the He li A4686 line 
strength. Therefore the assumption of BB central stars might introduce uncertainties in 
the determination of chemical abundances. We will also check in the next sections in a 
specific case the BB assumption vs. the use of a model atmosphere. We will see that 
differences in resulting abundances are quite small. 

The aim of our procedure is to match the intensity of the observed emission lines with 
that predicted by CLOUDY. The spirit of our procedure is similar to the one adopted by 
JC99 to determine the chemical abundances of PNe in M 31. In the following we briefly 
describe how our analysis works. 

CLOUDY needs the following input parameters: the central star energy distribution 
and luminosity, the nebular geometry, the electron density, and the chemical abundances. 
Because of the lack of information in our objects on the morphology and on the density 
distribution, we set the geometry to be spherical and the density to a constant value, 
although we will explore the effect of a r^^ density distribution in some PNe. The central 
star energy distribution was set to that of a BB. For the three PNe (PN91, 93, 101) with 
the best S/N spectra and with chemical abundances derived also with the ICFs method 
(M03), we will compare nebular abundances from the blackbody central star with those 
from the model atmospheres by Ranch l)2()()3|l (see Tabl^. 
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The first run is done setting the input parameters as follows. The BB tempera- 
ture is derived using the Ambartsumian's method f Ambartsumian I1932|l . i.e. using the 
He ii(4686)/H i(4861) line flux ratio, or the ([O iii](4959) + [0 iii](5007))/He ii(4686)) 
lines ratio method (Gurzadyan ^SS)) for medium and high excitation objects, respec- 
tively. The medium and high excitation is recognized from the observed line ratio 
He ii(4686)/H i(4861). Following Gurzadyan H199HI1 we caU high excitation PNe those 
with He ii(4686)/H i(4861) > 0.15. ^ 

The central star luminosity is set to reproduce the observed absolute flux of the 
[O III] A5007 nebular line, measured by Magrini et al. (|2U(J(J|I via aperture photometry 
and reddening corrected with E(B-V)=0.07 (van den Bergh H2000|l '). We made an effort to 
correct each A 5007 A [O in] flux with the individual reddening coming from the Balmer 
decrement. We saw that in some cases the corrected luminosity was substantially above 
the cutoff of the [O in] luminosity function. Thus we conclude that the de-reddening 
values from the Balmer decrement suffer from calibration issues. In these cases, while 
the line fluxes can be considered to be correct, c/3 is not appropriate to correct the 
photometric fluxes. We have then preferred to correct the photometric fluxes using the 
average extinction towards M 33. 

We varied the external radius from 10 ^^ cm (0.0032 pc) to 10 ^* cm (0.32 pc). This 
range corresponds to a reasonable interval of sizes of Galactic PNe. The electron density 
is evaluated from the [S n] AA6717/6731 line ratio. We have assumed the average non- 
Type I galactic PNe abundances from KB94 as initial chemical abundances and dust 
grains, with composition and size considered to be typical of Galactic PNe dust (cf. 
Fer land fl 99811 . The addition of dust in the model does not vary substantially the derived 
chemical abundances. 

Further iterations make the following adjustments. First the BB T^, is altered to match 
the He i/Hc n line ratio. When the He n A4686 line was too weak to be measured, we 
determined its upper limit so as to derive an upper limit to the central star temperature. 
These cases correspond to T* marked with b in Tabled while directly measured stellar 
temperatures are marked with a. In the few cases, where the [O n] A7325 doublet was 
measured, the nebular radius was varied to match the observed ratio of the nebular lines 
[O in] A5007 and [O ii] A 7325 (cases marked with c in Table 1) (cf. Che & Koppen lT!?^ 
In the other cases, models with different nebular radii were built. For these PNe, we found 
the smallest and the largest radius for which the model converged. Final quantities are 
the averages of quantities from models spanning over these nebular radius intervals (cases 
marked with d in Table 1). 



^ Gurzadyan (I1988II indeed distinguishes between low and high excitation objects. We prefer 
to use the term 'medium' in place of low because when the He ll A4686 line is measured, the 
excitation cannot be really regarded as low. We think this specification is more appropriate to 
PNe where the line is not detectable. 
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Subsequently chemical abundances are varied in order to match the observed and 
predicted intensities within 5% for the lines brighter than H/3 and 20% for the other 
lines. The Tg derived from the model is compared with the observed one, when available, 
i.e. in the three PNe where the [O in] 4363 A emission-line was measurable. In these 
PNe, the abundances of C and O (the most important coolants) were slightly modified 
to match the observed Te . The optical depth of the nebula is derived from the hydrogen 
ionization structure predicted by the CLOUDY model. The PNe are optically thick, while 
a few with lower density (PN28, PN96, PN125) resulted to be optically thin. 

Generally, we need from 10 to 30 iterations to reach convergence, i.e. to reproduce 
within 5-20% the observed line spectrum. Well determined quantities are the luminosity 
of the central star, because distance and absolute [O in] 5007 A flux are both well known, 
and the temperature of the central star (when both He i and He n lines are measured). 
The determination of the nebular radius is instead weaker, because the measurement of 
A 7325 A [O n] doublet is affected by large observational errors. In two PNe where the 
[O ii] A7325 doublet was measured (PN 91, 93) we introduced the dependence on r~^ in 
the density to find better agreement between observed and predicted [O iii]/[0 n] ratio. 
In fact the ratio [O in]/[0 n] is not only sensitive to the size of the nebula, but it depends 
also on several other parameters, mainly on the density distribution (cf. Che & Koppen 
I1983|l . When the central star temperature is derived from the upper limit of the He n 
A4686 line, only upper limit to the helium and other metal abundances can be derived. 

The derived central star temperatures and nebular quantities are shown in Tabled 
We recall that the T^, derived from the upper limit of the He ii A4686 are marked with b, 
whereas the directly measured T^, are marked with a. Derived nebular radii are marked 
with c, and upper limit nebular radii with d. 

In Table 12 chemical abundances for three PNe computed using ICFs (M03) are com- 
pared with current CLOUDY results. Note that for Ar/H the ionization correction factors 
are larger that one because only the [Ar in] A7135 line is detected. 

3.1. Tests of the procedure 

We have tested the accuracy of the method by applying it to seven bright Galactic 
PNe (Ml-8, M3-1, M3-5, NGC2867, NGC3195, NGC7009, NGC5852), using only the 
emission lines present in our extragalactic PN spectra, and comparing our results with 
the chemical abundances obtained with the ICFs method by KB94. PNe with central star 
BB temperatures spanning in a large range (from approximately 60,000 to 190,000 K) 
have been analyzed. For these tests we have used the basic approximations of our model: 
BB central star and constant density. The results are shown in Figures |2 and 13 We found 
a good agreement in the helium abundance, with a r.m.s. difference of 0.015 dex (triangles 
in FiglSJ. The typical r.m.s. differences of the two determinations, with CLOUDY and 
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with ICFs, of 0/H is 0.15 dex, of Ar/H is approximately 0.1 dex and of N/H is ~ 0.3 dex. 
In some cases these quantities are underestimated by our method whereas sometimes are 
overestimated. On the other hand S/H and N/H appear to be generally underestimated 
by CLOUDY, with a r.m.s. difference of 0.2 dex. A possible explanation for the systematic 
differences in S/H and N/H determination might be related to the fact that for these 
two elements only low ionization lines are observed. The considered PNe are medium 
and high excitation objects where these ions represent only a minor fraction of the total 
elemental abundance. Moreover the ionization potential of [S ii] is lower than that of 
[N ii], which might explain the larger discrepancy in S/H. 

4. Chemical Abundances 

We have applied the procedure above described to 11 PNc of M 33. In four cases we 
could determine the temperature of the central star within 5,000 K using He i and He ii 
lines. In the remaining cases we estimated an upper limit to the central star temperature 
using upper limits to the He ii A4686 flux. Considering the errors in the observed fluxes, 
the uncertainties in the model representation, and the quality in the flt obtained, we 
estimate typical errors are: in He/H of 0.05 dex (0.10 dex when T^, is an upper limit), 
of 0.2 (0.3) dex in 0/H, and of 0.3 (0.4) dex in Ar/H and S/H, and of 0.4 (0.5) dex in 
N/H. We note that S/H is generally under-estimated by the CLOUDY method in our 
test with Galactic PNe, and also N/H appears low compared to ICFs method for the 
highest values of N/H. 

Chemical abundances, central star parameters, nebular radius and electron density 
are shown in Tabled Chemical abundances previously obtained for the three PNe with 
the ICFs method (M03) are shown in Table [21 together with the new determinations, 
both with BB and model atmosphere f Ranch r2003|l . The four sets of determinations 
are consistent within the errors. Since the difference of the chemical abundances derived 
with CLOUDY using the two stellar ionizing sources are rather small, we used for the 
remaining PNe a BB central star. The use of a density dependent on r^^ with a BB 
central star produces better agreement between predicted and observed [O ni]/[0 ii] 
ratio. With this model we obtain a higher N/H and a lower 0/H, and consequently a 
higher N/0 than with ICFs method. 

S/H has been re-derived because of a numerical mistake in the previous work (M03). 
The discrepancy of S/H computed with ICFs and with CLOUDY (both BB and model 
atmosphere) is ~-0.2 dex, confirming the result mentioned above that CLOUDY tends to 
slightly underestimate S/H. 

The comparison in the brightest PNe of M 33 of abundances derived with the ICFs 
and those from photoionization models, confirms the accuracy of the CLOUDY procedure 
in measuring helium, oxygen, nitrogen and argon abundances. 
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Fig. 2. Comparison of He/H abundance of a sample of Galactic PNe derived with ICFs 
method by KB94 with He/H computed using CLOUDY. Triangles represent He/H when 
T^, has been derived in the present work. The id numbers refer to the Galactic PNe: 1- 
Ml-8, 2-M3-1, 3-M3-5, 4-NGC2867, 5-NGC3195, 6-NGC7009, 7-NGC5852. The derived 
Tj,/1000 is reported near to the corresponding measured He/H. 



Stellar evolution theory predicts that He and N are enhanced in PNe by nucleosynthe- 
sis process, whereas elements heavier than N, like O, Ar, S do not vary significantly from 
the moment of the formation of their central stars (cf. Iben & I},enzini llH83|l . In M 33, 
a comparison of our determinations of the chemical abundances of PNe (particularly for 
PNe whose T^ could be measured, see Table ^ with those of H ii regions (Vilchez et al 
I1988|l . especially for the best determined elements He/H and 0/H (Table QJ, confirms 
the reliability of the determinations. In fact, as expected, He/H (and also N/H, but re- 
member the larger errors that affect this determination) is enhanced compared to H ii 
regions while 0/H is almost unchanged. This indicates that He and N have undergone 
nucleosynthesis enrichment in PNe progenitors, as theory predicts. Consequently N/0 
appears to be increased in the M 33 PNe compared with H ii regions. 
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Fig. 3. Comparison of 0/H, N/H, S/H, and Ar/H abundances of seven Galactic PNe de- 
rived with the ICFs method by KB94 with the same quantities computed using CLOUDY. 

In Fig21 we present the derived oxygen abundances as function of the galactocentric 
projected distance. We find a possible consistency with oxygen gradient across the disc of 
M 33 detected by Vilchez et al. H1988|l for H ii regions. A weighted least squares fit gives 
a Alog(O/H)/Ad=-0.14 dex kpc~^, assuming the distance of 840 kpc, to be compared 
with the overaU gradient by Vilchez et al. lfTU55|l . Alog(O/H)/Ad=-0.10 dex kpc~\ 
considering the same distance. The Pearson correlation coefficient for our relationship is 
0.6. ON the other hand our result is quite dependent on a single PN (n. 8). Therefore 
the study of more PNe at high galactocentric distances is needed to clarify this aspect. 



4.1. Relationship among chemical abundances 

In Fig.|5jL we show the relation between He/H and N/0. Chemical abundances are plotted 
with their errors. The area included inside the dashed lines is the locus of Galactic PNe, 
derived from Fig. 3 in P04. Following the recent results by Exter et al. l2U(J4l for both bulge 
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Table 1. Central star temperatures and nebular parameters. Chemical abundances are 
expressed in 12+log(X/H). The line fluxes Fhq and i^[oin] are observed (MOO). The 
adopted reddening correction is E(B-V)=0.07 (van den Bergh l)2U(JU|l '). 
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(4) Dufour 119841 

(5) Leisy & Dennefeld i'}6tii\ 

(6) Dcnncfcld Jf{)ij{)l 



and disc Galactic PNe on N/O vs He/H, PNe with low N/O show also low He/H, but at 
N/O > 0.25 the whole range of He/H is sampled. The existence of a relationship between 
He/H and N/O has been discussed thoroughly by KB94, who found no correlation, 
although a correlation in the Type I PNe alone might exist. This issue is difficult to 
assess because there are intrinsic problems in determining helium abundance such as the 
collisional corrections and the uncertain ionization fractions. The data of M 33 for N/O 
vs He/H suggest a behaviour similar to the galactic one, but further studies are in order 
to really ascertain that. 

Fig-lit shows O/H vs N/O. The existence of an anticorrelation between these quan- 
tities is particularly controversial in the literature. Our data indicate a weak anti- 
correlation, with a Pearson correlation coefficient of 0.5 (0.2 in case of a weighted anti- 
correlation). Moreover, it appears that, for 12+log(0/H)<8.5, only high values of N/O 
are found, while there is a larger spread for higher O/H abundances. The existence of 
this effect has also been found in Galactic bipolar PNe (Perinotto & Corradi il998|l . 
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Fig. 4. The galactocentric trend of oxygen in M 33 from 11 PNe, recognized by their 
identification numbers in Tabled The solid line is the galactocentric gradient found from 
H II regions by Vilchez et al. H1988|l in the same galaxy. 



in the Galactic bulge (Stasinska et al. [TM5jl . in the SMC (Costa et al. IMM Leisy & 
Dennefeld lMiljI . possibly in M31 (JC99) and the LMC (Leisy & Dennefeld I^OTl) . al- 
though it is not apparent neither in the sample of Galactic PNe by KB94 (covering the 
range 8.2<12+log(0/H)<9.2) nor in bulge and disc Galactic PNe by Exter et al. I^UUl 
(8<12+log(0/H)<9.2). 

In the cases where this transition to high N/0 ratios for low oxygen abundances is 
well visible, it occurs at different 0/H for different galaxies, and roughly at the average 
overall metallicity computed from other classes of objects like HII regions. This suggests 
that the explanation offered by Costa et al. 1201)01 i.e. the dredge-up increasing N/0 in 
the AGB envelope is more efficient at low metallicities, might be not the main one, as the 
transition to high N/0 ratios would then occur always at the same metallicity. Together 
with the fact that in the Galaxy the effect is better visible for high mass progenitors like 
bipolar PNe (Perinotto & Corradi [l998|l . the data would instead point to a high efficiency 
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Table 2. Comparison among CLOUDY both with blackbody central star and a constant 
density, with model atmosphere by Ranch (|2003|l , and with BB and a density dependent 
on r^^ (two PNe) abundances determinations. The last raw gives the ICFs abundances 
determinations (M03). 



Id 


Method 


He/H 


N/H 


O/H 


S/H 


Ar/H 


91 


CLss 


11.31 


7.86 


8.61 


6.42 


6.26 




^i^rauch 


11.20 


8.01 


8.64 


6.52 


6.25 




CLss + dens-r-2 


11.31 


8.07 


8.48 


6.60 


6.18 




ICFm03 


11.28 


7.83 


8.70 


6.61t 


6.32 


93 


CLbb 


11.00 


7.53 


8.69 


6.46 


6.17 




^Ljrauch 


10.92 


7.55 


8.60 


6.44 


6.13 




CLss + dens-r-2 


10.93 


7.75 


8.47 


6.66 


6.09 




ICFm03 


11.00 


7.50 


8.85 


6.44t 


6.36 


101 


CLbb 


11.07 


7.33 


8.76 


6.18 


6.27 




CLirauch 


11.14 


7.34 


8.70 


6.13 


6.22 




ICFmo3 


11.15 


7.21 


8.72 


6.32t 


6.45 



fSee text. 

of the ON cycle in the most massive PN progenitors (Henry ^^nj, that would lower the 
O/H abundances while increasing N/0. 

In Figs.[5j: and Eli the relationship between N/0 and N/H are presented. In Fig. [S]; 
the dashed line represents the correlation found by Henry H1H9()|I for the Galactic PNe. 
Our data, shown in Fig. Et with their errors, show some scattering, but the Galactic 
relation between N/H and N/0 appears to apply as well to M 33 PNe. The Pearson 
correlation coefficient is 0.7. JC99 suggested the meaning of this correlation: its slope is 
close to unity and thus the variations of N/0 are mainly due to the variation of N/H 
(due to nucleosynthesis of N from progenitors of different mass), while oxygen abundance 
changes less during the life of PN central stars. In Fig. |3Ji, N/0 vs. N/H is reported for 
Galactic PNe (P04), for PNe of the LMC (Leisy & Dennefeld IMUjl and of M 33 (this 
paper). A similar trend is common to PNe of these three galaxies. 

We come now to the number of Type I PNe in our sample. The Type I PNe, according 
to the definition of Peimbert & Torres-Peimbert (|1983fl are those with log(N/O)>-0.3 and 
high He/H>0.125. A more stringent definition has been proposed by KB94, log(N/0)>- 
0.1. On the other hand it has to be noted that the definition of Type I PNe should be 
a function of metallicity, because the amount of N that can be produced by hot bottom 
burning is dependent on the amount of C present and also because the O/H abundance 
depends on the metallicity of the galaxy. Thus a different criterion should be adopted for 
Type I PNe in a galaxy of different metallicity from our Galaxy. The N/0 limit above 
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Fig. 5. Some relationships between chemical abundances, a) He/H vs. N/0. The area 
included inside the dashed lines is the main locus of Galactic PNe derived from Perinotto 
et al 120041 hereafter P04, Fig. 3. The dotted lines show the limits of type I PNe area, b) 
0/H vs. N/0: an anti-correlation in shown as seen by JC99 in M 31; M 33 PNe (present 
paper) appear to reach lower N/0 values than M 31 PNe. c) N/H vs. N/0. The dashed 
line shows the Galactic relation derived by Henry H1990fl as shown by JC99. In spite of 
some dispersion, the Galactic relation between N/H and N/0 appears to apply as well 
to M 33. d) N/H vs. N/0. Filled circles are PNe of M 33 (this work), triangles are PNe 
of LMC (Leisy et al. lM71| and squares are Galactic PNe (P04). 



which a PN is a Type I, meaning that it has experienced hot bottom burning conversion 
of C to N, should be the C+N/0 ratio from the H ii regions. This is difhcult to estimate 
since values for C are hard to obtain. Thus, since the metallicity of M 33 is not so different 
to that of our Galaxy, we still apply the Galactic criterion to discriminate Type I PNe. 
From Fig.Et, we note a lack in M33 of type I PNe, using the Peimbert & Torres-Peimbert 
l)1983|l criterion. Moreover, we find no Type I PNe with the definition by KB94. 
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The average log(N/0) of M 33 PNe is -0.89, see Fig.Eli. This value is lower than the 
Galactic non-Type I PNe (KB94) value (-0.55) and it is close to the non-Type I PNe 
value (-0.90; Leisy & Dennefeld l2004|l of the LMC, which is a galaxy of similar metallicity 
(see van den Bergh l^UUfl . whereas Galactic and LMC Type I PNe have a higher value 
(-f0.07 and -1-0.02, respectively, see Table [H). 

The number of Type I PNe might be connected to the metallicity of the host galaxy. 
In a deep survey for PNe of the SMC by Jacoby & De Marco H2002|l . they found an 
extremely large number of PNe where the [N ii] lines are stronger than Ha. This might 
translate to a large Type I PN fraction. Comparing the fraction of Type I PNe also in 
the Galaxy and in the LMC, they determined that there might be a trend of more Type 
I PNe for lower metallicity galaxy. They explained this trend with the fact that more 
massive PN central stars, which are progenitors of Type I PNe, produce more dust and 
thus end up being systematically dimmer and consequently harder to detect. This would 
be more so on high metallicity galaxy since more dust should be produced. M 33, which 
has a metallicity intermediate between the Galaxy and the LMC, might be expected to 
have a similar number of Type I PNe (between '^17% and ~20% cf. Jacoby & De Marco 

An explanation for the smaller number of Type I PNe detected might be that Type 
I PNe, which are principally represented by bipolar morphological type and generally 
associated with younger and more massive progenitors (cf. Corradi & Schwarz il995|l . 
are less luminous in the [O in] 5007 A line than other morphological types of PNe (e.g. 
Magrini et al. 120041 see Figs. 1 & 2, from which one notes that the cutoff of bipolar 
PNLF is fainter '^1.5 mag for LMC and ^2 mag for the Galaxy than that of elliptical 
PNe). We remind the reader that the PNe presented in this paper have been discovered 
with a survey complete up to ^2 mag fainter than the bright edge of the luminosity 
function (Magrini et al. I2000|l . If we consider that the results of luminosity functions of 
PNe with different morphology presented by Magrini et al. I2004L especially in the case 
of the Galaxy which has a similar morphology (M 33 is Sc III and the Galaxy is S(B)bc 
I-II, cf. van den Bergh l^OOfl . apply to M 33 as well, this means that a number of bipolar 
Type I PNe could have been lost because of the incompleteness of the survey; i.e. one 
should expose more in [O ill] to detect Type I PNe. 

4.2. Oxygen abundance and [O III] luminosity 

The determination of chemical abundances in an external galaxy rests on its brightest 
objects. We have examined the O/H vs [O in] luminosity in the M 33 sample, in the LMC, 
and in the Galaxy to test for possible biases in measuring O/H from the brightest PNe 
(see FigEJ. The sample of Milky Way PNe consists of objects in common between PNe 
whose chemical abundances have been re-determined by P04 and PNe whose distances 
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Fig. 6. [O III] luminosity vs 0/H for M 33 PNe (upper panel), MW PNe (central panel) 
and LMC PNe (lower panel). 

have been recently obtained by Phillips H2U(J2|I , whereas the LMC sample consists of PNe 
with chemical abundances and relative fluxes from the work of Leisy & Dennefeld H2004|l . 
and the absolute fluxes from Jacoby et al. p99(J|l . 

One could expect that there would be a correlation between the oxygen abundance 
and the [O in] 5007 A luminosity. However no trend of 0/H vs [O in] luminosity is seen in 
the three galaxies, as already noted for M 31 by JC99. There might be various reasons for 
that. The basic point is that oxygen is one of the most important coolants of the nebula 
so that a high O concentration lowers the nebular temperature and consequently the flux 
coming from the [O III] 5007 A line. There is therefore a feedback effect between a higher 
0/H concentration and a consequently lower Tg which produce the mentioned lack of 
trend of 0/H vs. [O in] luminosity. It means that 0/H derived from the brightest PNe 
(as from any PN) is representative of the total PNe population. This is particularly useful 
when deriving chemical abundances from PNe belonging to external galaxies, where only 
the brighter tip of their luminosity function is observed. 
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5. Summary and Conclusions 

Using the CLOUDY photoionization code, we have derived nebular parameters and some 
stellar effective temperatures for 11 PNe belonging to the spiral galaxy M 33, and mea- 
sured their He/H, 0/H, N/H, Ar/H and S/H abundances. The data are consistent with 
the 0/H gradient from H ii regions by Vilchez et al. H1988|l . but abundances in PNe far- 
ther from the centre of M 33 than the sample studied here are needed to ascertain a clear 
metallicity gradient from PNe in M 33. A trend of N/0 to N/H similar to the Galactic 
ones has been noted. We did not find PNe with the very high N/H and N/0 values that 
are typical of galactic Type I PNe. A possible explanation is that the incompleteness 
of the survey for M33 PNe (MOO) caused the loss of Type I PNe (N/O>-0.3), which 
are generally faint in the [O in] 5007 A emission line. Finally we found that the [O in] 
luminosity is clearly independent of the oxygen abundance. Therefore it is possible to 
use every PN, and particularly the brightest PNe, in external galaxies as representative 
of the whole PN population. 
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